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Knots and knotted fields enrich physical phenomena ranging from 
DNA and molecular chemistry to the vortices of fluid flows and tex- 
tures of ordered media. Liquid crystals provide an ideal setting for 
exploring such topological phenomena through control of their char- 
acteristic defects. The use of colloids in generating defects and 
knotted configurations in liquid crystals has been demonstrated for 
spherical and toroidal particles and shows promise for the devel- 
opment of novel photonic devices. Extending this existing work, 
we describe the full topological implications of colloids representing 
non-orientable surfaces and use it to construct torus knots and links 
of type (p,2) around multiply-twisted Mobius strips. 

knotted fields | topological defects | liquid crystals | colloids 

Controlling and designing complex three-dimensional tex- 
tures in ordered media is central to the development of 
advanced materials, photonic crystals, tunable devices or sen- 
sors and metamaterials [II El El HJ El El [71 [8j El [10] , as well as to 
furthering our basic understanding of mesophases [111 1121 [T3] . 
Topological concepts, in particular, have come to play an in- 
creasingly significant role in characterizing materials across 
a diverse range of topics from helicity in fluid flows [I4j [15] 
and transitions in soap films [16 , to molecular chemistry [17] . 
knots in DNA [18], defects in ordered media [19 , 20 , quantum 
computation [21, 22 and topological insulators [23]. Topolog- 
ical properties are robust, since they are protected against all 
continuous deformations, and yet flexible for the same reason, 
allowing for tunability without loss of functionality. 

Some of the most intricate and interesting textures in or- 
dered media involve knots. Originating with Lord Kelvin's cel- 
ebrated 'vortex atom' theory [24], the idea of encoding knot- 
ted structures in continuous fields has continued in magneto- 
hydrodynamics [25], fluid dynamics [15], high energy physics 
[26j [27l [28] , and electromagnetic fields [29 , 30 , and has seen 
recent experimental realizations in optics 31 , liquid crystals 
[32] and fluid vortices [33]. Tying knots in a continuous field 
involves a much greater level of complexity than in a necktie, 
or rope, or even a polymer or strand of DNA. In a field, the 
knot is surrounded by material that has to be precisely config- 
ured so as to be compatible with the knotted curve. However, 
this complexity brings its own benefits, for the full richness 
of the mathematical theory of knots is naturally expressed in 
terms of the properties of the knot complement - everything 
that is not the knot. In this sense, knotted fields are ideally 
suited to directly incorporate and experimentally realize the 
full scope of modern knot theory. 

Liquid crystals are orientationally ordered mesophases, 
whose unique blend of soft elasticity, optical activity and fluid 
nature offer a fertile setting for the development of novel meta- 
materials and the study of low dimensional topology in or- 
dered media. Much of the current focus centers on colloidal 
systems - colloidal particles dispersed in a liquid crystal host - 
which have a dual character. On the one hand, the liquid crys- 
tal mediates long-range elastic interactions between colloids, 
furnishing the mechanism for formation of colloidal structures 
and metamaterials [TJ [2] El HQ] • On the other hand, the col- 
loids, through anchoring conditions imposed by their surfaces, 
generate defects in the liquid crystal and so serve to induce 
and manipulate its topological properties. For instance, mul- 



tiple colloids exhibit a variety of entangled defect configura- 
tions [3H ES] , equally interesting states without defects [36] , 
and can even be manipulated so as to form arbitrary knots and 
links 32 , 37. . More recently, a significant advance has seen the 
fabrication of colloids with different topology [38] - tori up to 
genus five - verifying experimentally the relation between par- 
ticle topology and accompanying defect charge, and advancing 
a program to obtain topological control of materials through 
topological design. While the phenomena displayed by these 
systems is indeed rich, as surfaces these colloids (spheres, tori, 
etc.) all represent closed, orientable surfaces. 

In this article we extend these ideas to provide a complete 
topological characterization of all compact colloidal surfaces in 
a liquid crystal host. A/on-orientable surfaces fully exploit the 
non-orientable nature of liquid crystalline order [39] . We show 
that the topology of non-orientable surfaces enforces the cre- 
ation of topologically protected disclination lines. By varying 
the embeddings of the surface we exploit this topology to cre- 
ate metastable disclination loops in the shape of (p, 2) torus 
knots and links, for any p, around multiply- twisted Mobius 
bands. Through this combination of geometry and topology 
we elucidate a natural setting for the creation and control of 
complex knotted fields and the integration of mathematical 
knot theory into experimental science. 



Colloids, Surfaces and Topology 

Liquid crystals are typically composed of long, thin, rodlike 
molecules, which align in the nematic state along a common 
direction n called the director. This orientation is line-like 
rather than vectorial, meaning n ~ — n and the local ori- 
entation takes values in the real projective plane RP 2 , the 
ground state manifold for nematics [191 120] . The director 
field varies smoothly everywhere except on points or lines of 
discontinuity, known as topological defects, whose nature is 
captured by the way the order changes in their vicinity. Line 
defects, called disclinations, are characterized by the behavior 
on small loops around the singular line and so may be clas- 
sified by the fundamental group of the ground state manifold 
7Ti(RP 2 ) = Z2. Unique to liquid crystals, these defects reflect 
the non-orientability of nematic order with the molecules un- 
dergoing a 7r-rotation upon encircling the disclination. Simi- 
larly, point defects, referred to colloquially as hedgehogs, are 
classified by the second homotopy group 7T2(MP 2 ) = Z. This 
simple classification is augmented by the topological content 
of interactions between defects (characterized in terms of the 
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action of 7Ti on the higher homotopy groups, or by various 
Whitehead products [40]) so that, in particular, point defects 
of equal but opposite strength are equivalent (meaning freely 
homotopic) in the presence of a disclination and, whilst there 
is only one non-trivial element of tti(RP 2 ), disclinations that 
close up into loops fall into four distinct homotopy classes 
[4T] . Loosely, these may be thought of as corresponding to 
whether the linking number with other defects is even or odd 
and whether the loop carries an even or odd hedgehog charge 

Defects can be induced and manipulated by immersing col- 
loidal particles in the liquid crystal. This arises through the 
incompatibility of anchoring conditions on the particle sur- 
faces with the alignment imposed by the cell boundaries, or at 
large distances. The topological type of this incompatibility, 
or obstruction, can be associated with elements of the homo- 



topy groups 7T/c 



so that different surfaces can (loosely) 



be thought of as generating different types of defects. The na- 
ture of the obstruction depends on the anchoring conditions 
at the surface, but in the most common case of normal, or 
homeotropic, anchoring (as we consider here) it depends only 
on the topology of the colloid's surface. The classification 
theorem of surfaces [42 is a classic result of two-dimensional 
topology which states that any compact surface can be clas- 
sified up to homeomorphism by its genus, orient ability and 
number of boundary components. The genus is equal to the 
number of holes or handlebodies possessed by a surface; for 
example, a torus has genus one and a sphere genus zero. Ori- 
entability implies a consistent choice of normal vector can be 
made on a surface. The one-sided Mobius strip is the classic 
non-orient able surface; any normal vector on the strip will be 
flipped by going around the strip once, forbidding a consistent 
choice of surface normal. Finally, the number of boundary 
components is simply the number of distinct connected com- 
ponents in the surface boundary, e.g. a disk has one boundary 
component and a torus has none. 

While this is a complete topological classification of sur- 
faces in an abstract setting, for applications they must 
also be embedded (no self-intersections) into ordinary three- 
dimensional space, R 3 . Different embeddings are interesting 
in their own right - the whole of knot theory concerns embed- 
dings of a circle into R 3 - but they do not affect the homotopy 
class of the defect necessitated in the bulk, which we focus on 
first. Thus, with the classification of surfaces in mind it is 
natural to ask what topological implications each type of sur- 
face has for accompanying defects in the surrounding liquid 
crystal. The complete classification, summarized in Fig. 1, 
naturally separates into four classes of surfaces; orient able or 
non-orientable and closed or with boundary. 

Closed, orientable surfaces are known to induce defects 
corresponding to the element 1— g = x/2 of 7r2(RP 2 ), where 
g is the genus of the surface and x is the Euler characteris- 
tic [431 138] . Briefly, this relation comes through computing 
the degree of the Gauss map of the surface. The Gauss map, 
G, of a surface, X, is a map Q : X — » S 2 that sends every 
point of the surface to the direction of the surface normal at 
that point. For orientable surfaces with normal anchoring, 
the director can be given the orientation of the Gauss map, 
so that Q describes precisely the molecular orientation at the 
surface. The degree of this map - the number of times every 
point on S 2 is visited, counted with sign - is a homotopy in- 
variant [33] characterizing the type of defect that the surface 
generates 38 . Although experimentally the same surface can 
produce seemingly different defects, they are always charac- 
terized by this same element of tt2(RP 2 ). For instance, spher- 
ical colloids can nucleate either a point defect 1 or discli- 
nation loop [33] but the loop can always be shrunk continu- 
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Fig. 1. Topological characterization of compact surfaces with homeotropic bound- 
ary conditions, embedded in a three-dimensional, nematic liquid crystal. Closed ori- 
entable surfaces generate elements of 7T2(MP 2 ), equal to one minus the genus of 
the surface. Non-orientable surfaces with boundary generate the non-trivial element 
of 7Ti(RP 2 ), which forces the nucleation of disclination lines in the bulk. 



ously into a point [J5] , so that it is more properly classified by 
7T2(MP 2 ). More generally, orientable surfaces can never gen- 
erate elements of tti(RP 2 ), i.e. disclinations, as a topological 
requirement. Their orient ability ensures that any disclination 
loops formed can always be removed in pairs or shrunk into 
points. 

Closed, non-orientable surfaces cannot be embedded in R 3 
without self-intersection, meaning that a true representation 
of any of these (e.g. the real projective plane or the Klein 
bottle) in a liquid crystal is not possible. For this reason, we 
do not consider closed, non-orientable surfaces any further. 

Orientable surfaces with boundaries have trivial topologi- 
cal implications for the surrounding liquid crystal. Their ori- 
ent ability, as in the closed case, forbids them from generating 
elements of tti(RP 2 ) as a topological necessity. In addition, 
since they have a boundary, they cannot generate any ele- 
ments of tt 2 (RP 2 ). We argue as follows. Closed, orientable 
surfaces separate space into an inside and an outsid^] and 
they generate point defects in both regions. If one cuts a hole 
in the surface, creating a boundary component, then these 
defects can be combined to leave a defect-free texture. Cut- 
ting more holes in the surface will not change things as one 
could always make the texture on the new hole identical to 
the surface normal that was removed. 

Non-orientable surfaces with boundary necessarily gener- 
ate a non-trivial element of tti(RP 2 ). A generalization of the 
Gauss map to non-orientable surfaces Q : X — >• RP 2 assigns 
to every point of the surface the line element (point in RP 2 ) 
corresponding to the direction of the (unoriented) surface nor- 
maQ Now consider the loop space of our surface, QX. Com- 
position with the non-orientable Gauss map Q creates a set of 
representatives of tti(RP 2 ). If X is non-orientable then there 
must be at least one map in this set which represents the 
non-trivial element (and generates a disclination). Suppose 
there was no such map, then every map in Q [QX] could be 
lifted from RP 2 to S 2 and we would have created an orientable 
Gauss map, Q, implying the surface is orientable, a contradic- 



*We consider only "tame" embeddings such that the piecewise linear Schonflies theorem applies, 
t It is more common, and more useful in general, to think of the non-orientable version of the Gauss 
map as a map to the Grassmannian Gr2(K 3 ) - see, e.g. Milnor JW, Stasheff JD (1974) Char- 
acteristic Classes (Princeton Univeristy Press, Princeton), pp 55-70 - but since Gr2(R 3 ) is 
canonically homeomorphic to MP 2 there is no loss in our discussion. 
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tion. The disclinations created in this way must entangle the 
surface, since any disk spanning a non-orient able loop on the 
surface must be pierced by the defect. 

With closed, orient able surfaces, different elements of 
7T2(MP 2 ) could be generated according to the genus of the sur- 
face. In the same way, we might ask if different non-orient able 
surfaces with boundary can act to generate different kinds of 
disclination loops in the bulk liquid crystal. Although there 
is only one non-trivial element of 7Ti(RP 2 ), there are four dis- 
tinct homotopy classes of disclination loops [411 120] (the four 
types may be thought of as corresponding to even/odd link- 
ing number with other disclinations and even/odd hedgehog 
charge). We must, therefore, determine which of these is gen- 
erated by the surface. We cannot force the existence of linked 
loops in a path connected domain. If we have a connected 
domain then a single disclination loop can go through all the 
non-trivial cycles created by the surface. Since we will then 
have only one disclination, it cannot be linked. Furthermore, 
it must have zero hedgehog charge. The surface simply acts 
to smoothly align the director field in some region of space so 
that the texture is equivalent to one containing a lone discli- 
nation loop, without the colloid but with the same director 
orientation in its place, which must have zero charge. Thus all 
non-orientable surfaces with boundary, independent of their 
topological type, have the same topological implication for the 
liquid crystal, in contrast to the closed, orientable case. These 
four classifications are summarized in Fig. 1. 



Mobius Strips and Knotted Defects 

Non-orientability of the surface enforces the existence of a 
disclination loop but it leaves open the precise form of the 
defects and their equilibrium configuration. These are deter- 
mined by energetics and by the nature of the embedding of 
the surface. While surfaces with boundary are unavoidably 
two-dimensional, experimentally realizable, but topologically 
equivalent, surfaces may be constructed by using thin mate- 
rial with homeotropic boundary conditions on the faces and 
planar anchoring on the thin edges, as shown schematically 
in Fig. 2. In this way a colloid with varying surface an- 
choring conditions can be made to faithfully represent a two- 
dimensional surface. Such boundary conditions are essential 
to mimic non-orientable surfaces and ensure the topological 
properties we describe - fully homeotropic boundary condi- 
tions simply replicate an orientable torus. Modern fabrication 
techniques allow for the manufacture of such exotic surfaces 

EHJ. 

As the Mobius strip is the prototypical non-orientable sur- 
face - all non-orientable surfaces contain the Mobius strip as 
a subset - it serves as an elementary guide to the behav- 
ior of non-orientable surfaces in liquid crystals. A Mobius 
strip with homeotropic boundary conditions will generate a 
non-trivial element of 7Ti(RP 2 ) as one passes around the strip 
and hence must be threaded by a disclination loop of zero 
hedgehog charge, entangling the surface. The shape of discli- 
nation loops around colloidal particles is governed largely by 
the requirement to minimize the distortion energy of the sur- 
rounding director field. A simple heuristic for this can be 
constructed as follows. As one passes through a disclination 
line a rotation of approximately tt/2 is induced in the director. 
This rotation mediates the transition from the surface normal 
orientation to that of the far-field and will best minimize the 
distortion in the director field when it is concentrated along 
those parts of the surface where the local anchoring and far- 
field directions are perpendicular. Since we consider colloids 
with homeotropic anchoring this simply gives the requirement 




Fig. 2. Simulation results of a Mobius strip with homeotropic boundary conditions 
on the flat faces, shown schematically in black along the strip. Note the planar an- 
choring on the short edge required to correctly represent a two-dimensional surface. 
The two-dimensional surface can then be thought of as living in the center of the 
strip. The defect line (red) is clearly visible entangling the strip, as predicted by JT|. 
The inset shows the local director profile near the defect, showing the +1/2 twist 
disclination on the inside of the strip; there is a corresponding -1/2 twist profile on 
the outside. 

that 

S n -n o = 0, [1] 

along the disclination, where S n is the surface normal and no 
is the far- field director. For a spherical particle in a uniform 
far-field, this predicts that the disclination will lie on a great 
circle in a plane perpendicular to this far-field direction, which 
is the observed position of Saturn ring defects [33]. Likewise 
the twisted shape of disclinations around spherical colloids in 
a cholesteric [3SJ [37] is correctly predicted by the same heuris- 
tic. 

The preferred defect configuration for a Mobius strip can 
be found by numerical simulation using continuum Landau-de 
Gennes modeling (see Methods). As shown in Fig. 2, the min- 
imum energy configuration is a single disclination loop entan- 
gling the strip, in the location predicted by Q. Of course, the 
precise configuration depends on the strip's orientation rela- 
tive to the far- field, but we have found that the orientation 
shown, with the strip's centerline in a plane perpendicular to 
the far- field orientation, has the lowest observed energy. The 
cross-section of the disclination loop shows a twisted —1/2 
profile on the outside of the strip and a +1/2 twisted profile 
on the inside, as it has to in order to carry no hedgehog charge 
[47] . This hedgehog charge may be computed by several meth- 
ods. The recently developed Pont ryagin- Thorn construction 
[13] and methods related to disclination profile switching [33 
both assert that the charge of the disclination is zero, as re- 
quired on topological grounds. 

Perhaps the simplest generalization of the Mobius strip 
topology is to vary its embedding in R 3 . Different embed- 
dings can be obtained by changing the number, p, of half- 
twists that the strip contains. The 'canonical' Mobius strip 
has one half- twist (p = 1), but more generally if p is odd then 
the surface is still non-orientable and has the same topology as 
the Mobius strip. When p is even the surface is orientable and 
has the topology of an annulus. Nonetheless, the embeddings 
are distinct and carry their own topological embellishments. 
The boundary of a Mobius strip is a circle. For a single half- 
twist this is a simple unknot, but for p half- twists it is a (p, 2) 
torus knot (p odd) or link (p even). To see this, note that the 
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Fig. 3. Knotted and linked disclinations in chiral nematics stabilized by the pres- 
ence of twisted surfaces with homeotropic boundary conditions on the flat faces, and 
planar on the thin edges. They are torus knots and links of the form (p, 2). Shown 
are: (a) p = 2 Hopf link, (b) p = 3 trefoil knot, (c) p = 4 Solomon's knot and 
(d) p = 5 cinquefoil knot. The defects in the center are hedgehogs, existing in pairs 
above and below the strip. 

boundary of a p-twisted strip (p odd) lives on a torus whose 
major radius is that of the strip and whose minor radius is half 
the strip width. The curve the boundary draws on this torus 
goes round the meridional cycle p times - once for each half- 
twist - while traversing the longitudinal cycle twice, which is 
the definition of a (p, 2) torus knot. The story is the same 
for orient able strips with p even, except that there are two 
components to the boundary and they form a link. 

Can this structure, coming from the nature of the embed- 
ding, be exploited to controllably produce knotted and linked 
disclination loops in liquid crystals? Disclination lines that 
follow the surface of such a multiply-twisted Mobius strip will 
have the same shape and properties as the colloid boundary, 
yielding precisely constructed knots and links. Here we show 
that such configurations can be stabilized in chiral nematics; 
examples for doubly, triply, quadruply and quintuply twisted 
strips are shown in Fig. 3. They produce the p — 2 Hopf link, 
the p — 3 trefoil knot, the p — 4 Solomon's knot and the p — 5 
cinquefoil knot, respectively. All these knots and links obey 
the topological requirements set out in the first section; the 
strips with an even number of twists are topologically trivial 
and the strip is entangled by an even number of disclinations, 
those with an odd number of twists are non-orientable and so 
enclose an odd number of disclinations. 

The stablization of these knotted structures is not just a 
question of topology, energetics also come into play. In this 
regard, the chirality (inverse pitch) of the system has an im- 
portant role in the stability of these configurations. Indeed 
it is generally true that chiral systems allow for more exotic 
structures [371 14]. In the achiral nematic system, the knotted 
defects are unstable and the liquid crystal assumes a ground 
state configuration consisting of p small disclination loops en- 
tangling the strip along the contours where S n - no = 0, as pre- 
dicted by JTJ. A similar configuration, with slightly twisted 
loops (Fig. 4), is also the ground state in cholesterics - the 




Fig. 4. Comparison of configurations for a quintuply twisted strip, (a) Ground- 
state configuration consisting of 5 small disclinations loops at the locations predicted 
by [fljl They have a slightly twisted profile, caused by both the twisting of the strip 
and the chiral nature of the system, (b) (5, 2) torus knot entangling the strip, with 
two hyperbolic hedgehogs nucleated above and below. The energy difference between 
these configurations is approximately 1.2%. 

knots are metastable - although the difference in energies is 
small (of order 1-2%) and decreases both with increasing chi- 
rality and knot complexity p. The behavior with increasing p 
can be understood in terms of the total length of disclination 
line, which scales as p for the isolated loops and as y4 + p 2 
for the knots. If the chirality is increased such that the pitch 
becomes smaller than the width of the colloid then the discli- 
nations develop twists analogous to those around spherical 
colloids [461 E3- 

Like cholesterics, torus knots have a handedness - the 
(p, 2) and (p, —2) knots are mirror pairs - and it is not surpris- 
ing to find that the relative handedness influences the stabil- 
ity of the textures. The handedness of the knot is set by the 
shape of the colloid, and if this matches that of the cholesteric 
then knotted textures are stable, otherwise the system tries 
to expel the reverse twist in the configuration, resulting in an 
unstable knot. 

The disclination lines themselves represent just a small 
portion of the entire system. Their knottedness imprints a 
complex orientational order on everything that is not the knot. 
Fig. 5a shows a slice through the director field of a quintu- 
ply twisted Mobius strip in a cell with fixed normal boundary 
conditions. To match the boundary conditions of the cell, the 
knotted disclination is accompanied by two hyperbolic hedge- 
hogs, expanded into small loops, above and below the strip. 
The disclination-colloid pair has a constant profile that simply 
rotates uniformly as one moves around the strip; qualititavely, 
this structure does not depend on p, the number of half-twists. 
The cross-section through the colloid has a profile reminiscent 
of a double twist cylinder [25], split apart by the colloid into 
two separate +1/2 twist disclinations. 

These configurations bear a striking similarity to the re- 
cently discovered 'toron' textures [4 . Torons come in several 
flavors, but the ones that concern us have two hyperbolic de- 
fects above and below a double twist torus; that is, a double 
twist cylinder - the building block of blue phases [35] - bent 
around such that it forms a torus. Like the knots of Fig. 3, 
these textures are stabilized in thin cells of chiral nematics, 
with uniform alignment at large distances. More than a su- 
perficial similarity, the knotted field configurations we present 
here can be thought of as the result of expanding the central 
part of a toron locally into two +1/2 twisted defects, that 
then rotate to form knots. The whole configuration is then 
stabilized by the colloid; indeed if the colloid is artificially 
removed from the simulation after equilibration, the config- 
uration collapses into a toron. This connection may provide 



* Since the hedgehogs are actually small loops in the simulations, an equivalent procedure is to take 
those loops as additional boundaries of a Seifert surface for the link consisting of the torus knot and 
two unlinked unknots. The genus of this Seifert surface is the same as that of the other one. 
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Fig. 5. Field Structure, (a): Director field profile of a (5, 2) torus knot and the two accompanying hyperbolic hedgehogs, opened up into small rings, above and below the 
strip, with disclinations bound to the colloid's surface. The local disclination-colloid profile is largely constant up to a rotation as one moves around the strip. Though this 
profile is for a p = 5 configuration, qualitatively the director field does not change as one varies p. The far-field director is fixed to be vertical on the boundary of the cell, 
(b): Pontryagin-Thom surface [13] for the trefoil knot-colloid pair. The surface is constructed by considering all points where the director field is perpendicular to the far-field. 
The surface is then colored with the remaining RP 1 degree of freedom (as illustrated in the inset). The defects, disclination line and two accompanying hedgehogs, are shown 
in black. The color winding around the hedgehogs establishes them as unit charge defects. The surface can be patched by removing the hedgehogs, and allowing the surface 
to pass through the colloid. One then obtains a Seifert surface for the knot, composed of two disks (top and bottom) connected by 3 (generally p) twisted bands. It is readily 
verified that this surface has genus 1 (generally (p — l)/2, for p odd), the genus of the trefoil knot. While there is also a color winding around the disclination, this is a 
four-fold winding, giving the disclination even (equivalent to zero) hedgehog charge, as required. 



a potential route for the production of these textures using 
Laguerre- Gaussian polarizing beams. 

A more quantitative investigation of the texture can be 
made through use of the Pontryagin-Thom (PT) construc- 
tion [13], allowing us to examine the full topology of the field 
structure and even perform rudimentary 'experimental knot 
theory'; we can use the texture around the knot to compute 
the knot genus, a classical knot invariant 49 . The PT con- 
struction involves drawing the surfaces on which the director 
field is perpendicular to a given direction: such a surface is 
shown in Fig. 5b. Surfaces constructed this way must end 
either on topological defects or, as in our case, colloids. For a 
knotted disclination, this PT construction then yields a Seifert 
surface for the knot, and it is a classic result that the minimal 
genus over all possible Seifert surfaces is a knot invariant for 
the given knot [50 . To obtain a proper Seifert surface, the PT 
surface must be continued through the colloid and the points 
corresponding to the hedgehogs filled irj^] Since our disclina- 
tions are torus knots, they have genus (p — l)/2 for p odd and 
(p — 2)/2 for p even [50]. While, in principle, the genus of 
the PT surfaces is only bounded from below by this number, 
in the simulations presented here the limit is always reached, 
and we may read off the correct genus. 

In conclusion we have demonstrated the design and con- 
struction of elegant knotted and linked disclinations in liq- 
uid crystals, exploiting natural embeddings of twisted strips 
and a complete classification of the topological implications 
of homeotropic colloids. These knotted disclinations impose a 
complex structure on the director field, allowing knot invari- 
ants to be computed from the texture and fascilitating the 
full introduction of knot theory into experimental soft matter 
systems. The configurations themselves are intimately related 
to the recently discovered toron textures, suggesting a route 
for experimental realization. Future work will further explore 



the value of knot theory in understanding the subtleties of 
knotted fields and their potential for novel devices. 

Materials and Methods 

To simulate the twisted strip colloids we use the standard Landau-de Gennes Q-tensor 
formalism. The liquid crystal order parameter is taken to be a traceless, symmetric 
tensor field Q(r). The free energy of the system, F , is then given as 

+ |tr((V x Q + 2<7o<2) 2 ) jd 3 r 

+ 1™ tr((Q-Q°) 2 )dV. 

[2] 

The domain Q, is the volume of the liquid crystal, and S is the surface of the colloid. 
The constant Aq defines the bulk energy scale of the material and 7 represents 
an effective temperature. L, the elastic constant, controls elastic distortions in the 
director field. ^0 is the chirality of the liquid crystal and is equal to 2tv / P where 
P is the cholesteric pitch. ^0 > denotes a right-handed material and the pitch 
was set to be of order the colloid diameter, though the configurations were found to 
be stable at higher chiralities. Finally, W controls the strength of the anchoring on 
the colloid surface, with the preferred orientation given by Q°. The free energy is 
minimized by evolving the Q-tensor according to Landau-Ginzburg relaxational dy- 
namics. These equations are then solved using finite difference methods on a cubic 
mesh of size 250 X 250 X 250 or 260 X 260 X 150. Simulation parameters 
(A 0.16, L 0.03, 7 3.2, W 0.2) were chosen to match typi- 
cal values for a cholesteric liquid crystal with pitch comparable to the colloid radius 
(varied between 1 and 2 times this) and strong surface anchoring. 
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